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We study a generalization of the voter model on complex networks, focusing on the scaling of 
mean exit time. Previous work has defined the voter model in terms of an initially chosen node 
and a randomly chosen neighbor, which makes it difficult to disentangle the effects of the stochastic 
process itself relative to the network structure. We introduce a process with two steps, one that 
selects a pair of interacting nodes and one that determines the direction of interaction as a function 
of the degrees of the two nodes and a parameter a which sets the likelihood of the higher degree 
node giving its state. Traditional voter model behavior can be recovered within the model. We 
find that on a complete bipartite network, the traditional voter model is the fastest process. On a 
random network with power law degree distribution, we observe two regimes. For modest values of 
Q, exit time is dominated by diffusive drift of the system state, but as the high nodes become more 
influential, the exit time becomes becomes dominated by frustration effects. For certain selection 
processes and parameters values, an intermediate regime occurs where exit occurs after exponential 
mixing. 

PACS numbers: 89.75.Fb, 02.50.Ey, 89.75.Hc 

I. INTRODUCTION 

The voter model has been extensively studied on lattices [1] and, in recent years, on complex networks [2-6] and 
is closely related to models of language evolution [7], ecological dynamics [8], opinion dynamics [9], and epidemic 
spread [10]. The voter model defines a dynamical process where nodes arc each assigned one of two states, +1 or 
— 1. Connections are defined on a lattice by nearest neighbors or on an arbitrary network by edges. Each update step 
consists of selecting a pair of nodes and giving the state of one node to the other. In the most frequently studied 
version of the model, the first node chosen adopts the state of the second. The most interesting object of study is 
the mean exit time, i.e. the mean time to achieve complete agreement. For the voter model on complex networks, 
a node will be chosen second with frequency proportional to its degree, and so its influence is fixed by the selection 
process. In this paper, we introduce a generalized voter model with a single tunable parameter that allows control of 
the influence of topology in a manner independent of the selection process. In this generalized model, the probability 
of a node giving its state to its neighbor is proportional to k", where k is the node's degree and a is chosen. 

Voter model processes based entirely on selection frequency are denoted either link update or node update [4]. In 
link update dynamics, every time step an edge is selected uniformly at random. One of the two nodes at the ends of 
the link is then chosen randomly to give its state to the other. In node update dynamics (described above) a node 
is selected to adopt the state of a random neighbor. The neighbor chosen at random is likely to have high degree, so 
high degree nodes have more influence. 

II. GENERALIZED VOTER MODEL 

The two processes described above specify both the selection of a pair of interacting nodes and which node adopts 
the state of the other. We separate this process into two distinct steps to better understand the contribution of each 
aspect of the process. 

Given a network of size N, each node i has state Sj = ±1 and degree fc^. We deflne Pjj to be the probability of 
giving to node i the state of node j during a given time step. There are two independent components of this event: 
the probability Sij of selecting an edge connecting nodes i and j and the probability W{ki, kj) that a node with degree 
kj gives its state to a node with degree fcj. Thus, 

P,,=S,,Wih,k,). (1) 

The form of W{ki, kj) is motivated by comparison with node update dynamics in uncorrelated networks. In uncorre- 
lated networks, node update dynamics can be described by considering all nodes of like degree to be indistinguishable 

in the ensemble average [3]. In the following, {«} and {7} refer to the subgroups of all nodes with degrees ki and kj, 
respectively. Under the node update process, the probability of giving a node with degree k^ the state of a node with 
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degree kj is 



Pij=n,'^ (2) 



where n,, nj are the fraction of nodes with degrees fcj and kj and Hi = J2i ''^i^i average degree. This can be 

interpreted as the probabiUty rii of selecting a node with degree ki times the probabiUty of following an edge into a 
node with degree kj. A node in {j} has kj edges, so the probability of following an edge into it is proportional to kj. 
If node pairs are selected in this manner, a particular pair of nodes in {i} and {j} can be chosen by either picking 
from {i} and following an edge to a node in {j} or by picking from {j} and following an edge to a node in {i}. For 
node selection, this gives 

S^j='^{k, + kj). (3) 

We propose a generalization that includes the standard node update dynamics, which requires SijW{ki,kj) = Pij. 
The only form of W{ki, kj) to do this is 

This form is also consistent with our definition of W{ki,kj) as a probability. Notably, it suggests a one parameter 
generalization; 

W{ki,kj,a)=j^^. (5) 

Qualitatively, the parameter a determines how much a node asserts its degree when transmitting its state. For 
a > 0, the higher degree node of a pair is more likely to give its state to the lower degree node, a bias that increases 
with a. For a < 0, the opposite is true. The special case a = ignores topology in determining the direction 
of interaction since W{ki,kj,0) = 1/2 always. Edge update dynamics is recovered by using edge selection to find 
pairs and setting a = 0. Node update dynamics occurs when node selection determines pairs and a = 1. A recently 
investigated "invasion" dynamic, where a node is picked to give its state to a random neighbor (opposite the traditional 
model), occurs for node selection and a = —1. If all nodes of the network have the same degree, as in a mean field or 
lattice topology, then all values of a are equivalent to the traditional voter model. 

Effectively, this model assumes some connection between the behavior of the agents and the underlying network on 
which they live. For example, if this were to be thought of as a model of opinion dynamics, a value of a > 1 under 
node update selection would correspond to a situation where an individual prefers to behave like those who are more 
connected than himself. Celebrities or well-regarded experts are extremely infiuential, for example, but the same 
forces that drive their high visibility also keep them from being influenced by the non-famous. The forces of influence 
and accessibility compete, such that a small value of a makes all nodes able to change state quickly, but limits the 
influence of any node. A high value of a makes high degree nodes influential, thus able to order their neighborhood 
quickly, but those influential nodes will flip only on rare occasions. 



A. Dynamics 

To understand the dynamics, we study the master equation for an arbitrary network. The probability of a system 
being in state s — {si} at time t is defined to be P{s,t). Denote by the state s where Si i—^ —Si and let Sij be 
the probability of selecting the edge between nodes i and j. For brevity, we write W{ki,kj,a) = Wij. The master 
equation is 

|p(s, i) = S,,W,j (^^) P(s\ t) - S,jW,j (^^) P(s, t) (6) 

Let pk be the ensemble average probability of a node with degree k being in a +1 state. Averaged over all random 
graphs, nodes of like degree are indistinguishable and we consider each such like-degree subgraph separately [3]. The 
evolution of arbitrary ensemble average functions can be found in a straight-forward manner [11], giving: 



(7) 
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We can find a conserved magnetization, p*, by choosing coefficients Cj such that 

fc 

Since Sij is symmetric in i and j, this can happen for arbitrary subgraph densities only if 



C.^ = G^. (9) 
rife ni 



This implies that Cfe cx nfe/c". Normalizing, 



P = 



= — J^nfcfcVfe (10) 



where is the a*'' moment of the degree distribution. Note that the ensemble conserved magnetism is independent 
of the process of selecting node pairs. 



III. BIPARTITE NETWORK 



The simplest degree homogeneous topology is the fully connected bipartite network. Such a network is given by two 
groups of nodes, group A with size a and group B with size b. A node in group A is connected to every node in group 
B, but none in group A. The degree of nodes in A, Ua, is the size of B, giving kA = b and similarly ks = a- In this 
situation all edges are interchangeable, so there is no difference between the two selection processes. We will consider 
only the effect of a. Let pa be the concentration of +1 opinions in A and pt, be the concentration of +1 opinions in 
B. In our model, the special value a = 1 is equivalent to the case studied in [3] on the same network and we follow a 
similar procedure, omitting details that can be found there. Prom Equation (10), the conserved magnetization is 



= ^.-1 + ^.^1 (fe-Va + a-V) (11) 

For any initial conditions, the ensemble average subgraph densities approaches p*. If all nodes in A start as +1 and 
all nodes in B start as —1, then the probability of ending in the +1 state is 

P+ = 6a-i+aa-i (12) 
The mean exit time is given by the backward Komologorov equation [12]. solves 

-^^^ = (Pa - PbW-^db - a^-^da)T+ ^{pa + Pb ~ 2paPb){a^-^dl + b^~^d^,)T. (13) 

where da and dh are partial derivatives with respect to the initial subgraph densities. The first term describes convection, 
which brings the subgraph densities to some equal value, and the second term describes the diffusion of the network- 
wide state [3]. The convective dynamics can be shown to be fast for all a. 

The fast step toward equal subgraph densities has a negligible impact on extinction time and we can consider only 
the subsequent one dimensional problem. We define p = Pa = Pb and apply a change of variables using Equation (11). 
After integrating, 

T = -{a'- + 6i-«)(a«-i + b^-')^{plog{p) + (1 - p) log(l - p)). (14) 

a + b 

This has a similar form to the standard voter model, but with a factor that is symmetric about a = 1 and non- 
vanishing. If we take a = Xb, then 

Toe (2 + Ai-" + A"-i)-^6. (15) 
1 + A 

If A ^ 1, corresponding to a star-like graph, 

r~Al"-^l6. (16) 

This scaling is confirmed in simulations (see Figure 1). Notably, the standard voter model, a = 1, is the fastest 
process for any complete bipartite network. 
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IV. ARBITRARY NETWORKS 



Similar analysis extends naturally to networks in which a node's degree determines the network structure. Many 

random network models fall into this category, notably any random network generated by the configuration model, 
including those with scale-free distributions, and Erdos-Renyi networks [13]. Small world networks are not included, 
however, as certain nodes have exceptional topological characteristics that are independent of their degree [14] . 

As in Equation (13), we can write the equation satisfied by the mean exit time on an arbitrary network. In the 
following. Pi refers to the density of +1 states on the subgraph of nodes with degree ki = i. 

St = Y^ SijWtjiPj - Pi)SAT +IY1 SijWtjiPi + Pj - 2pipj)6fdfT (17) 

ij ij 

The system is again split into a convective term and a diffusive term. The assumption of fast approach to well-mixed 
state must be treated more carefully in our generalized model, but there do exist cases where diffusion dominates. 
Namely, it was observed in [3] that node update dynamics on a scale free network has fast convection compared to 
its diffusive exit time. There is, in general, a prefactor to the term pi{l — pj) when evaluating the probability that 
two nodes have differing states [6] . Numerical simulations show that this is of order one and not generally a constant 
over time when a ^ 1, so it does not affect the final scaling results. 

We expect this to not be true for all a. In the case of a » 1, a node with degree higher than all its neighbors 
will act to dictate its neighbors' states, but only rarely be changed itself. The network in this case may not be able 
to quickly approach the global equilibrium given by p* , since these locally highest degree nodes will be pinned for a 
time dependent on a. If this duration is longer than the time for the rest of the system to become ordered via mixing 
and diffusion, a quasi-frustrated state occurs. The exact local topology, rather than just degree distributions, can 
dominate the dynamics. 

Let us suppose that the exit time is diffusion dominated and will return to discuss the validity of this assumption. 
The system can be approximated by a one dimensional equation in p = ^ rijfcf p^: 



-N : 

And thus 



^ {y. j f (1 - pHT. (18) 

"-EWF- '''' 

The denominator can be simplified by noting that 

Si^w^^kr = \[Z ^nW^j^i'^ + ^ Si^w^k^'^ ) = ^ ? 

ij \ ij ij J ij 

This gives: 

l^ij '^ij'^i i^j 

Since Sij is a probability, the sum can be thought of as a weighted average over selection probabilities. Interestingly, 
for a = 0, neither the form of interaction selection nor the network topology matter. In that case, po = 1 and 
Wij = 1/2, so 

T„=o « N. (22) 

This agrees with the observation in [2] that exit times scale with N in situations where the unweighted magnetization 
is conserved, which corresponds exactly with a = 0. 

To go farther, we need to specify the selection scheme and the network. We focus our consideration on random 
uncorrelated scale-free networks with degree distribution Uk ~ k~'^ . Networks with power law distributions appear in 
a variety of social and biological contexts and exhibit a range of interesting behaviors [13] . Let us first consider node 
update, for which Sij = riirij ^2^^^ . Then Equation (21) becomes: 

TvcxiV^. (23) 
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The a*'' moment can be approximated by an integral: 



k°'n{k)dk (24) 
up to an effective maximum degree kmax-, defined by n{k)dk = 1/N [15]. It is easily 

{N^ a> v-1 
N^""^^ u-2 <a <u-l (25) 
N a<u-2 

For V > 2 and any a, the exit time increases without bound as system size increases. We simulated the process 
on random network generated by the configuration model [16] and found good agreement with our predictions (see 
Figure 2). 

For edge update dynamics, Sij = nitij ^1^. _ Lqw degree nodes are selected less frequently under edge selection than 
node selection. The diffusive exit time can be calculated similarly, giving: 



Te=N('^] . (26) 



The approximate scaling for edge update is 




a> v — 1 

T£<x{ n'''^ v-2<a<v-l (27) 
a <v -2 

This leads to very different scaling behavior. For the parameter regions 

i/<3 a>v-l (28) 

^ < ^"'"'"^ u-2<a<v-l (29) 
o 

the diffusive exit time vanishes as N increases. Simulations show that there is a diffusive region which agrees with 
our predictions for smaller values of a (see Figure 3). 

The convective process involves an exponential decay of each pi to its equilibrium value, with rate determined by 
the network structure. When this value is comparable to l/riiN (or 1 — l/riiN) for all i, the system effectively reaches 
convergence before diffusive time scales matter. We can read off the dynamics from Equation 7. The decay rate r is 
not strongly dependent on system size, so the relationship with A'' comes only from the exit condition 

1 ~ e-^^'\ (30) 

This implies that when diffusive time scales vanish, 

r£~log(iV). (31) 

Exponential mixing is only observed as a dominant time scale under link update selection for values of A'', v, and a 
such that frustration does not occur. 

To understand where frustration comes from, consider a node whose degree is higher than that of all of its neighbors. 
In the limit a ^ oo, such a node will take an arbitrarily long time to flip. For large values of a, escape from this 
frustrated state must then be the dominant time scale, as diffusive time scales for high a are dependent only on v. 
For cither selection process, the probability of a given node being sampled in a time step is a constant, independent 
of system size. It is also possible for small connected subgraphs to be frustrated if every edge out of the subgraph 
is to a node of lower degree. As a result of locally highest degree nodes and clusters, a small number of hard-to-flip 
nodes can hold back a large network from reaching convergence (Figure 4). 

The continuum treatment in Equation 17 averages over network ensembles before solving for the exit time. In the 
high a case, considering an averaged network ignores local leader effects and fails to give an accurate solution. The 
approach to an ensemble equilibrium state does not occur, so diffusion about this state is not a valid assumption 
(Figure 5). 
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The degree distribution of locally highest degree nodes can be approximated for the non-assortative case quite 
simply. The degree distribution of local leaders, pu{k), is the independent product that a node has degree k and that 
all k neighbors have a degree less than k [17]: 

Mk)=p{k)(Y.'^] . (32) 
\k'<k ^1 / 

No term is strongly related to system size, so the total number of local leaders Nu scales linearly with N. 

The dynamics of these nodes are based on extremely local behavior and thus very hard to approximate, but we can 
construct the slowest possible node and understand its behavior. For a sufficiently large system, one such node will 
likely exist and thus dominate the exit time. The probability on any given time step to flip a local leader with degree 
I is formally given by: 

Pfiipil) = E ^i^k^Pi^i + (33) 

i * 

where P{(Ji ^ ui) is the probability that the neighboring state differs from the local leader's state and subscripts 
index individual nodes. 

For either selection process, when I is large, 

Sii-\. (34) 

We assume the probability of a neighbor being in the opposite state is finite and treat the following as a lower bound. 
Simulations show that this is reasonable. Assuming that I" ^ k", 



m^^^\^y.ki\- (35) 

The slowest situation occurs when the local leader in question is surrounded by relatively low degree nodes and 
P{1) l^~"/N. For this case, the time to flip in units of system size scales as: 

Tf ~ r-\ (36) 

The highest degree of the local leaders scales linearly with the global highest degree [17], so for < 3 we have 
Tf ^ ]\[('^~^)/(''~^) , This will be true for large a, but requires the system size to be high enough that at least one rare 
node such as the one described is likely. Such networks are prohibitively large to sample, but the qualitative situation 
described explains why frustration sets in at a smaller value of a fov v = 2.4 than u = 2.8 in Figure 2. Local leaders 
for small v have a higher degree than for large v, though numerical comparison from Equation 32 shows them to be 
rarer. 



V. SUMMARY 



Recent work [18] has found an approximate mean exit time for a duplication process on networks with arbitrary 
edge weights, assuming that diffusion is the dominant time scale. In this work, we demonstrate that there are at least 

two natural ways for this estimate of exit time to fail. As observed in Baxter et al., the time for the system to reach 
a metastable equilibrium can be at least as large as the diffusive exit time scale. We see this in the edge selection 
process for values of a and v where the diffusive exit time vanishes as system size gets large. The frustrated dynamics 

in the node selection mode, however, presents a new way in which the diffusive estimate can fail. System dynamics 
are driven by a small number of topologically special nodes, breaking the assumption that a continuum description 
applies. 

We have defined and analyzed a single parameter voter-model-like stochastic process that is identical to the original 
voter model on a homogeneous network. On a complex network, our generalized voter model has a tunable dependence 
on local network topology, allowing us to control the differing effects of selection and the direction of influence. On 
complete bipartite graphs, the traditional voter model is the fastest process to reach an absorbing state. On scale 
free networks, the dynamics depend strongly on the selection process. Node selection has two regimes; a diffusive 
one characterized by a well deflned average magnetism and diffusion constant based on global network properties, 
and one with frustrated dynamics stemming from the neighborhood around a small number of locally highest degree 
nodes. Edge selection, previously considered uninteresting, has three regimes. In addition to diffusive and frustrated 
regimes, it also has a middle ground characterized by exponential mixing. Understanding the dynamics involved in 
this transition to frustration would be an interesting avenue for future work. 
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FIG. 1: Circles are the simulated fit of dlogT / d\og N for a complete bipartite network with M = 40 and A*" ranging from 100 
to 5000. The dotted line is the scaling of Equation 16 for the same range of N. 
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FIG. 2: Simulated values of dlogT / dlog N under node selection based on several hundred runs for A'^ from 750-15000. Circles 
correspond to v — 2.8, squares to = 2.4. The dashed line is the scaling based on the diffusive estimate, calculated by fitting 
the numerically calculated sum in Equation 21 for similar values of A'^. Note that frustration begins to dominate for u = 2.4 at 
a > 1.6, causing the deviation from the diffusive estimate. 
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FIG. 3: Simulated values of rflogT/dlog A'^ under edge selection from several hundred runs of values of A'^ from 750-15000. 
Circles are for i/ = 2.8, squares for f = 2.4. The dashed line comes from fitting the numerically calculated sum in Equation 21 
for similar values of A''. The horizontal line is the effective slope of T ~ log A*' for the system sizes used. 
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FIG. 4: A example of node update dynamics on a small network containing a global highest degree node (G) and a separate 
local highest degree node (L) for a = 1 and a = 10. The darkness of a node corresponds to the average fraction of time spent 
in a state opposite the final state of the network over 1000 realizations with identical initial conditions. A darker node has 
spent more time in a contrary state than a light node. For a = 1, states are well mixed and T = 9.4. For a = 10, mixing does 
not occur. The local leader and its neighborhood spend most of the time in a contrary state and T = 4814. 
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FIG. 5: Typical dynamics of the fraction of active edges (edges connecting nodes with different states) for diffusion and 
frustration. Both have N = 10000, v = 2.8, and edge selection. The top line, in black, is for a = 1 and the bottom line, in 

gray, is for a = 8. The a = 1 case is diffusive and fluctuates to convergence after reaching the ensemble average value. The 
a = 8 case decays exponentially to a value greater than the ensemble average value, because a locally highest degree node or 
cluster is slow to flip. Just before t — 1200, these nodes flip and the system reaches the absorbing state. Inset shows early time 
dynamics on a similar network for a = 10 in gray and a = 4 in black. Both decay exponentially initially, but the high a case 
becomes frustrated and the other continues to convergence. 



